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Die Erfindung betrifft ein Verfahren zur Erfassung dynamischer 
Systeme, die durch im Zeitverlauf nichtstationare 
Systemparameter charakterisierbar sind, insbesondere ein 
Verfahren zur Segmentierung von Zeitreihen von Mefigrofien 
(Variablen) dynamischer Systeme und zur Identif izierung der 
die Segmente charakterisierenden Systemparameter (Moden) . 

Unter einem dynamischen System wird hier insbesondere jede 
Erscheinung verstanden, deren Zeitverlauf durch eine diskrete 
Abbildung vom Typ 



x(t+U = f am (a(t)) (O.i) 

darstellbar ist. Es werden jedoch auch Systeme mit mehreren 
(z. B. zwei) simultan auf genommenen Zeitreihen x, y gemaB 

y(t+x) = f a(t) (x(t)) (0.2) 

betrachtet. Dabei bezeichnen a(t) einen Satz 
charakteristischer Systemparameter, x einen Zustand, der im 
allgemeinen einen Vektor in einem mehrdimensionalen 
Zustandsraum bildet, und y einen zeitlich verschobenen 
Zustand. Der Zustandsraum wird von Variablen aufgespannt, die 
z. B. physikalische, chemische, biologische, medizinische, 
geologische, geometrische, numerische und/oder 
prozeBtechnische GrolJen sein konnen. 



Die Anzahl der Systemvariablen, die zusammen mit der Dynamik f 
das System beschreiben, entspricht der Dimension des 



Zustandsraumes. Hier werden Systeme betrachtet, deren 
Parameter a ebenfalls zeitlich veranderlich sein konnen. Ein 
gegebenes System mit zeitlich unveranderlichen Parametern a 
wird im folgenden auch als Mode bezeichnet. 

Beobachtbare oder mefibare Systemvariablen (MeBgroflen) bilden 
erfaUbare Zeitreihen oder Datenstrome, die fur die jeweilige 
Abfolge von Systemmoden charakteristisch sind. Sind die 
Systemparameter innerhalb der Zeitreihen liber bestimmte 
Zeitabschnitte (Segmente) unveranderlich, so kann die 
Zeitreihe entsprechend den jeweils bestehenden Systemmoden 
unterteilt werden (Segmentierung) und jedes Segment einer 
Systemmode zugeordnet werden (Identif izierung) . 

Viele Vorgange in der Natur wie auch bei technischen 
Anwendungen wurden sich dann vorhersagen und/oder 
kontrollieren lassen, wenn die ihnen zugrundeliegenden 
dynamischen Prozesse mathematisch modelliert werden konnen. 
Die Analyse und Charakterisierung praktisch gegebener 
dynamischer Systeme wird oft dadurch erschwert, daii sich die 
Systemmoden wahrend der Beobachtung andern. Beispiele hierfur 
sind allmahliche Veranderungen, die sich in Drifts oder Trends 
der Systemparameter widerspiegeln, oder auch spontane oder 
abrupte Veranderungen in der Dynamik komplexer Systeme, etwa 
wenn sich Konf igurationen spontan oder von auJien getrieben 
plotzlich andern. 

Ein Beispiel fur ein betrachtetes System ist die Erzeugung von 
Sprachsignalen im Mund-Rachen-Raum, bei der das System standig 
seine Konf iguration und somit seine Mode andert. Es besteht 
ein starkes Interesse daran, die Moden zu erfassen und zu 
identifizieren, die einem beobachteten Zeitverlauf einer 
Variablen (im Beispiel: Luf tdruckschwankungen) zugrundeliegen, 
auch urn verbesserte Vorhersagen iiber das betrachtete System 
treffen oder dieses besser kontrollieren zu konnen. 



Prinzipiell konnen dynamische Systeme anhand vom gemessenen 
Signalen analysiert werden und es sind Reihe von Methoden 
bekannt, aus Zeitreihen Modelle zu gewinnen, die sich fur die 
Vorhersage und eine Kontrolle des Systemverhaltens eignen. So 
ist bekannt, dafi sich der Zustand eines dynamischen Systems 
durch Erfassung der Zeitabhangigkeit beobachtbarer Mefigrofien 
modellieren laflt. Diese Modellierung erfolgt gemaiJ einem 
ersten Ansatz durch eine Rekonstruktion des Zustandsraumes 
mittels sogenannter Zeit-Verzogerungs-Koordinaten, wie es z. 
B. von N. H. Packard et al. in "Physical Review Letters" {Bd. 
45, 1980, S. 712 ff.) beschrieben ist. Auf der Basis einer 
solchen Rekonstruktion lafit sich dann nur ein einziges 
(globales) Modell f fur die Dynamik finden. Die globale 
Rekonstruktion des Systems ist ferner nachteilig, da bei 
Anwendungen fur vieldimensionale Systeme eine Vielzahl von 
EingangsgroBen als Randbedingungen vorab bekannt sein mussen 
und/oder aufgrund der hohen Dimensionalitat das System sich 
praktisch nicht mehr schatzen (erfassen, abbilden) lafit 
und/oder ein ubermafiig hoher, nicht praktikabler Rechenaufwand 
entsteht . 

Aufierdem ist dieses Verfahren fur den Fall zeitlich 
veranderlicher Parameter allgemein nicht anwendbar. Die 
Analyse und Modellierung von dynamischen Signalen wird jedoch 
haufig dadurch erschwert, dafi die zugrundeliegenden Systeme 
sich in der Zeit in wesentlichen Parametern verandern. 
Beispiele sind z. B. Signale aus der Medizin, bei denen ein 
Organ wie das Herz oder das Gehirn viele dynamische Moden hat, 
die sich abwechseln, oder etwa den Sprachsignalen, bei denen 
das erzeugende System, der Mund-Rachen-Raum, in der Zeit 
offensichtlich verschiedene Konf igurationen durchlauft. 

Gemafi einem weiteren Ansatz ist daher aus der Publikation von 
K. Pawelzik, J. Kohlmorgen und K.-R. Muller in "Neural 
Computation" (Bd. 8, 1996, s. 340 ff.) bekannt, Datenstrome 
entsprechend zunachst unbekannten, zeitlich wechselnden 
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Systemmoden durch Simulierung mit mehreren, untereinander in 
Konkurrenz stehenden Modellen zu segmentieren. Die Modelle 
werden bevorzugt durch neuronale Netzwerke gebildet, die 
jeweils fur eine Dynamik charakteristisch sind und nach 
vorbestimmten Trainingsregeln urn die Beschreibung der 
einzelnen Punkte des Datenstroms konkurrieren. 

Mit dieser Methode gelingt es, eine Zeitreihe in Abschnitte 
quasistationarer Dynamik zu zerlegen und simultan Modelle fur 
diese System-Moden aus den Zeitreihen zu identif izieren. 

Die Segmentierung gemaii K. Pawelzik et al., zu der unten 
Einzelheiten angegeben werden, erlaubt die Zuordnung von 
Segmenten zu bestimmten Systemdynamiken oder -moden und ftthrt 
zu einer Erfassung des Datenstroms als Vorgang mit einem 
diskreten Schalten (sog. "switching") zwischen den Moden. 
Diese Beschreibung der Parameterdynamik komplexer Systeme 
stellt zwar gegenuber der o. a. globalen Modellierung 
hinsichtlich der Genauigkeit und der Segmentierung 
verschiedener Systemzustande einen Fortschritt dar. Allerdings 
kann der Ubergang zwischen verschiedenen Systemzustanden nicht 
ausreichend beschrieben werden. Es hat sich insbesondere bei 
der Analyse realer Systeme z. B. bei medizinischen Anwendungen 
gezeigt, dafi die Segmentierung auf bestimmte Falle mit 
moglichst klaren Modenunterschieden und geringem Rauschen 
beschrankt, allgemein jedoch bei zeitlichen Veranderungen der 
erzeugenden Systeme unzuverlassig ist. 

Solche zeitlichen Veranderungen der erzeugenden Systeme machen 
die beobachtbaren Signale nichtstationar und fiihren dazu, dafi 
sich die Systeme im allgemeinen nicht mehr durch einheitliche 
Modelle beschreiben lassen. Erfolgen diese Veranderungen der 
Systeme plotzlich, so spricht man von Sprungprozessen. 

Die Aufgabe der Erfindung ist es, verbesserte Verfahren zur 
Erfassung der Moden dynamischer Systeme mit nichtstationaren 
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Systemparametern anzugeben, mit denen die Beschrankungen 
herkommlicher Verfahren iiberwunden werden konnen und die es 
insbesondere ermoglichen, automatisch mit praktikablem 
Bearbeitungsaufwand und hoher Zuverlassigkeit Segmentierungen 
und Identifizierungen von Zeitreihen mit einer erhohten Zahl 
von Einzelheiten vorzunehmen. 

Diese Aufgabe wird durch das Verfahren mit den Merkmalen gemafi 
den Patentanspruch 1 gelost. Vorteilhafte Ausfuhrungsformen 
der Erfindung ergeben sich aus den abhangigen Ansprtichen. 

Die Erfindung basiert auf der Idee, Obergange zwischen 
verschiedenen Moden eines dynamischen Systems als 
Zwischenmoden des Systems aufzufassen, die paarweise lineare 
Interpolationen der Ausgangs- bzw. Endmoden des Obergangs 
darstellen. Die betrachteten dynamischen Systeme gehen eher 
allmahlich von einer Mode in die andere iiber, statt abrupt 
zwischen Moden zu schalten. Die Erfindung zielt darauf ab, 
solche Obergange zwischen dynamischen Moden in Signalen und 
die Moden zu identif izieren. 

Es wird daher bei einem Verfahren zur Erfassung der Moden 
dynamischer Systeme z. B. nach einer Schaltsegmentierung einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
Systems eine Drif tsegmentierung vorgenommen, bei der in jedem 
Zeitabschnitt, in dem das System von einer ersten Systemmode s L 
zu einer zweiten Systemmode ubergeht, eine Folge von 
gemischten Vorhersagemodellen g^^ erfaiit wird, die durch eine 
lineare, paarweise Oberlagerung der Vorhersagemodelle f t j der 
zwei Systemmoden s l#j gegeben ist. 

Gegenstand der Erfindung ist auch eine Vorrichtung zur 
Erfassung eines dynamischen Systems mit einer Vielzahl von 
Moden Si mit jeweils charakteristischen Systemparametern a(t) . 
Die Vorrichtung enthalt eine Einrichtung zur Aufnahme einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 



Systems, eine Schaltsegmentierungseinrichtung, die dazu 
eingerichtet ist, in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fur die Systemvariablen x(t) ein vorbestimmtes 
Vorhersagemodell f t fur eine entsprechende Systemmode Si zu 
erfassen, und eine Drif tsegmentierungseinrichtung, mit der in 
jedem Zeitabschnitt, in dem das System von einer ersten 
Systemmode s ± zu einer zweiten Systemmode s 3 ubergeht, eine 
Folge von gemischten Vorhersagemodellen gi erfaflt wird. Die 
erfindungsgemafie Vorrichtung kann ferner Einrichtungen zur 
Einstellung von Interpolations- und Segmentierungsparametern, 
Vergleicherschaltungen zur Verarbeitung von Vorhersagef ehlern 
von Vorhersagemodellen, Anzeige- und Signalisierungsein- 
richtungen und Speichereinrichtungen enthalten. Die 
erfindungsgemafle Vorrichtung kann ein Monitor fur 
physiologische Daten oder physikalische oder chemische 
ProzeJiparameter sein. 

Mit der Erfindung steht ein Instrument zur Verfugung, welches 
ein grofies Anwendungspotential in vielen medizinischen, 
wissenschaftlichen und technischen Bereichen hat. Mit der 
Segmentierung von Signalen und einer entsprechenden 
Identifizierung der zugrundeliegenden Dynamik eroffnen sich 
neue Moglichkeiten der Vorhersage und der Kontrolle auch 
wesentlich nichtstationarer Systeme. 

Anwendungen der Erfindung haben gezeigt, dafi sich 
kontinuierliche Ubergange zwischen Systemmoden sicher 
identifizieren lieBen und daii die zugrundeliegenden Dynamiken 
durch die Modelle mit einer Prazision beschreibbar sind, die 
in vielen Fallen eine Vorhersage des Systemverhaltens 
ermoglicht. Mit der Erfindung lassen sich in vielen Fallen 
nichtstationarer Prozesse Modelle identifizieren, die sich 
auch fur eine Kontrolle der Prozesse eignen, die ohne eine 
Berilcksichtigung der Nichtstationaritat nicht moglich waren. 
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Ausfiihrungsformen und weitere Vorteile der Erfindung werden i 
folgenden unter Bezug auf die beigeftigten Zeichnungen 
beschrieben. Es zeigen: 

Fig. 1: Kurvendarstellungen zur Illustration eines ersten 
Segmentierungsschrittes des erf indungsgemafien Verfahrens; 

Fig. 2: Kurvendarstellungen zur Illustration eines weiteren 
Segmentierungsschrittes des erf indungsgemafien Verfahrens 

Fig. 3: Kurvendarstellungen zur Segmentierung von 
Blutregulierungsdaten nach dem erf indungsgemafien Verfahren; 
und 



Fig. 4: Kurvendarstellungen zur Segmentierung von EEG-Daten 
nach dem erf indungsgemafien Verfahren. 

Im folgenden werden zunachst Einzelheiten der Erfindung unter 
Bezug auf die Figuren 1 und 2 und anschliefiend praktische 
Anwendungsbeispiele erlautert. Es ist dem Fachmann 
ersichtlich, dafi die Erfindung nicht auf die 

Anwendungsbeispiele beschrankt ist, sondern entsprechend auch 
in anderen Gebieten angewendet werden kann, wie sie 
beispielhaft weiter unten aufgefuhrt sind. 

(1) Erfassung von Drif tubergangen in nichtstationaren 
Zeitreihen 



Erfindungsgemafi werden nichtstationare Zeitreihen mit einer 
zweischrittigen Prozedur erfafit, bei der zunachst eine 
geeignete Modellierung und anschliefiend eine sogenannte 
Drift segmentierung erfolgt. Die Modellierung ist dazu 
eingerichtet, in jedem Zeitabschnitt einer vorbestimmten 
Mindestiange fur jeden Systemparameter ein vorbestimmtes 
Vorhersagemodell fiir eine entsprechende Systemmode zu 
erfassen. Hierzu erfolgt vorzugsweise eine herkommliche 
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Schaltsegmentierung, wie sie beispielsweise aus der 
Publikation von K. Pawelzik et al. in "Neural Computation" 
(Bd. 8, 1996, s. 340 ff.) bekannt ist. Die Modellierung kann 
aber auch durch eine andere, in Bezug auf die gewonnene 
Systeminformation zur Schaltsegmentierung aquivalente Prozedur 
erfolgen, die an eine konkrete Anwendung z. B. bei bekannten 
reinen Moden oder Randbedingungen angepaJit ist. 

Die Schritte der Schalt- und Drif tsegmentierung werden im 
folgenden naher erlautert. Hinsichtlich der 

Schaltsegmentierung wird hiermit der Inhalt der Publikation 
von K. Pawelzik et al. ausdriicklich vollstandig in die 
vorliegende Beschreibung einbezogen. 

(i) Schritt 1 (Schaltsegmentierung) 

Die Schaltsegmentierung dient der Ermittlung 
charakteristischer Prediktoren, die zur Beschreibung der 
Systemmoden geeignet sind. Die Schaltsegmentierung kann 
entweder an einer Trainings-Zeitreihe oder an der zu 
untersuchenden Zeitreihe durchgefuhrt werden. In beiden Fallen 
konnen die ermittelten Vorhersagemodelle oder Prediktoren fur 
weitere, unbekannte Zeitreihen verwendet werden. 

Es wird ein dynamisches System mit einer endlichen Zahl N 
unterschiedlicher Moden betrachtet. Fur die j-te Mode ist ein 
Wert (allg.: Vektor oder Satz) a,(t) eines beobachtbaren 
Systemparameters charakteristisch, der jeweils mit einer 
Funktion f i(t) (i = 1,...,N) aus einem Satz von N Funktionen f 
modelliert werden soil. Es wird die Zeitreihe {x t } = Xj (t) der 
Systemvariablen betrachtet und zeitabhangig jeweils nach der 
Funktion f i(t) gesucht, fiir die {y t } = Yj (t) = fi^fx^t)) eine 
neue Zeitreihe vorherzusagender Punkte y 3 (t) darstellt, die in 
Bezug auf die Systemmoden qualitativ die selben Eigenschaf ten 
hat wie {x t }. Durch den Wechsel der Modellfunktion f im 
Zeitverlauf wird die Schaltsegmentierung gefunden, die die 
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Zeitreihe {x t } entsprechend der wechselnden Systemmoden 
unterteilt . 

Die Funktionen f werden als Prediktoren (oder: 
Vorhersagemodelle, Expertenfunktionen) aus einem Satz von 
Netzwerken mit veranderlichen Parametern durch ein geeignetes 
Trainingsprogramm abgeleitet, bei dem simultan sowohl die 
Parameter der Netzwerke als auch die Segment ie rung ermittelt 
werden. Der Begriff "Netzwerk" wird hier fur alle moglichen 
geeigneten Model If unktionen verwendet, also vorzugsweise fur 
neuronale Netzwerke, aber auch z ., B . fur Polynome oder lineare 
Funktionsapproximationen. Die optimale Wahl eines neuronalen 
Netzwerkes erfolgt in Abhangigkeit von der spezifischen 
Anwendung. Es werden vorzugsweise Netzwerke mit einer 
schnellen Lernf ahigkeit, wie z. B. sogenannte RBF-Netzwerke 
(Radial Basis Function Network) vom Moody-Darken-Typ 
verwendet . 

Das Training erfolgt unter der Voraussetzung, dafi die 
Systemmoden nicht mit jedem Zeitschritt wechseln, sondern eine 
geringere Schaltrate aufweisen, so daS eine Systemmode tlber 
mehrere Zeitschritte erhalten bleibt. Die angenommene Grenze 
der Schaltrate bzw. Anzahl der Zeitschritte, uber die eine 
Systemmode erhalten bleibt, ist zunachst ein freier 
Eingangsparameter und kann anwendungsabhangig in geeigneter 
Weise beispielsweise in Abhangigkeit vorgegebener 
Erfahrungswerte oder nach einer Parameteranpassungsstrategie 
ausgewahlt werden. Bei der Parameteranpassungsstrategie kann 
vorgesehen sein, fur die Schaltrate einen Ausgahgswert 
vorzugeben und mit diesem einen Vorhersagefehler (siehe unten) 
zu ermitteln. Falls die Schaltrate zu hoch oder zu niedrig 
gewahlt ist, fiihrt eine entsprechende Uber- oder 
Unterspezialisierung zu einem zu hohen Vorhersagefehler. Im 
weiteren Verlauf der Anpassung kann dann die Schaltrate 
optimiert werden, bis der mittlere Vorhersagefehler unterhalb 
vorbestimmter Grenzen liegt. 
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Das Training erfolgt durch Maximierung der Wahrseheinlichkeit 
W, dafi der Satz von Netzwerken die Zeitreihe {x t } erzeugen 
wiirde. Es handelt sich urn ein Training mit konkurrierendem 
Lernen (engl.: "competitive learning"), wie es im einzelnen in 
der Abhandlung "Introduction to the theory of neural 
computation" von J. Hertz et al . (Addison-Wesley Publishing 
Company 1991, insb. Kap. 9: "Unsupervised competitive 
learning") beschrieben wird. Aus dieser Abhandlung ist die 
anwendungsabhangige Implementierung eines solchen Trainings 
ableitbar. Die Trainingsregel des konkurrierenden Lernens auf 
der Grundlage des beim Lernen auftretenden Fehlers ist gemafi 
(1) darstellbar. 



Slog W 

"af— * 



e -P(y-fj) 2 

2 je -P<y- f i) 2 



(y - f i> (1) 



Diese Trainingsregel stellt sicher, daJi die Lerngeschwin- 
digkeit (Verbesserung der Parameter) fur die Funktionen f mit 
geringstem Abstand zum Zielwert y am hochsten ist. 

Figur 1 zeigt das Ergebnis der Schaltsegmentierung am Beispiel 
der Analyse einer chaotischen Zeitreihe {x t } mit x ttl = f (x t ) 
zwischen den vier Moden: 

f x (x) = 4x(l-x) fiir x e [0, 1] 

f 2 (x) = f^f^x)) 

f 3 (x) = 2x fur x <= [0, 0.5) bzw. 

f 3 (x) = 2(l-x) fur x g [0.5, 1] 

f 4 (x) = f 3 (f 3 (x)) 



FUr die ersten 50 Zeitschritte wird mit einem Startwert 
x 0 = 0.5289 zunachst f : angewendet. Anschliefiend erfolgt ein 
Ubergang (Einzelheiten siehe (ii) ) zur Mode f 2 , die nach 
Schritt 100 bis Schritt 150 stationar wird. Entsprechend wird 
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ab Schritt 200 bzw. Schritt 300 jeweils fiir 50 Schritte die 
Mode f 3 bzw. f 4 eingenommen. Danach erfolgt der Rticktibergang 
auf f x . Figur la zeigt einen Ausschnitt (Schritte 300 bis 450) 
des Zeitverlaufs der Zeitreihe {x t } mit x t+1 = f(x t ). 

Die Segmentierung der ersten 450 Zeitschritte mit 6 
Prediktoren f^, i = 1,...,6 (RBF-Netzwerke vom 
Moody-Darken-Typ) ist in Figur lb dargestellt. Das Training 
ergibt eine Spezialisierung von vier der Prediktoren (6, 2, 4, 
3) jeweils auf eine der o. a. vier Moden. Die stationaren 
Bereiche liegen bei den Intervallen [0, 50] und [400, 450] 
(f x ), [100, 150] (f 2 ), [200, 250] (f 3 ) und [300, 350] (f 4 ) . Die 
ubrigen zwei Prediktoren (3, 5) haben sich auf die 
Obergangsbereiche zwischen den Moden spezialisiert . Dies zeigt 
den Nachteil der herkommlichen Schaltsegmentierung, bei der im 
Fall von Obergangen der entsprechende Zeitbereich ohne 
adaquate Beschreibung mehrfach unterteilt wird. 

Anstelle des hier beschriebenen sog. "Hard competition"- 
Trainings, bei dem bei einem Trainingsschritt jeweils nur ein 
Vorhersagemodell optimiert wird ("winner takes all") kann auch 
vorgesehen sein, den Grad der Konkurrenz im Rahmen eines sog. 
"Soft competition"-Trainings zu verandern, wie es im einzelnen 
in der Publikation von K. Pawelzik et al. beschrieben ist. 

(ii) Schritt 2 (Drif tsegmentierung) 

Beim zweiten Schritt werden die Obergange (sog. Driften, 
nicht-abruptes Obergehen, gleitendes Wechseln) zwischen den 
Systemmoden berucksichtigt . Mit der Erfindung wurde als 
wichtige Voraussetzung fur die Drif tsegmentierung 
herausgefunden, dafi der Obergang von einer ersten Systemmode 
direkt in eine zweite Systemmode und nicht Uber eine dritte 
Systemmode erfolgt. Das Driften zwischen Systemmdden wird 
somit wie folgt durch eine Oberlagerung von (oder paarweise 
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lineare Interpolation zwischen) genau zwei Moden modelliert . 
Dabei treten gemischte, ggf. abgestufte Zwischenmoden auf, die 
jedoch nicht^eigene (reine) Systemmoden sind. 

Es wird ein Satz von P reinen Systemmoden, die jeweils durch 
ein Netzwerk k(s) , s e P, reprasentiert werden, und ein Satz 
von M gemischten Systemmoden betrachtet, die jeweils durch 
eine lineare Uberlagerung von zwei Netzwerken i(s) und j (s) , 
s sM, reprasentiert werden. Das Modellnetzwerk g 3 fur eine 

gegebene Mode seS, S = PuM wird gemaB Gleichung (2) 
dargestellt. 

_ j f k(sj <x t ) fur s g P 

g s (x t ) ^\^(s)f i(s) (x t ) + b(s) f j(s) (x" t ) fur sg M (2) 

In (2) ist it der Vektor (x t/ x t _ T , . . . , x tHn _ 1Jt ) der 
Zeitverzogerungskoordinaten der Zeitreihe {x t } und f ifj sind 
Prediktoren, die gemaft der o. a. Schaltsegmentierung ermittelt 
wurden. Dabei ist m eine Einbettungsdimension und x der 
Verzogerungsparameter der Einbettung. Die Einbettungsdimension 
ist die Dimension des Phasenraumes, in dem das System 
betrachtet wird und in dem die Modelle operieren. 

Fur jede gemischte Systemmode sind zwei Parameter a, b 
zusammen mit zwei Netzwerkindizes i, j charakteristisch. Zur 
Vereinfachung des Rechenaufwandes wird die Zahl der gemischten 
Moden begrenzt. Mit 0 < a(s) < 1 und b(s) =1 - a(s) wird 
eine endliche Zahl von Werten a(s) definiert. Zur weiteren 
Vereinfachung werden gleiche Abstande zwischen den Werten a(s) 
gemaii Gleichung (3) gewahlt. 

a r = r / (R +1) mit r = 1,..., R (3) 

R entspricht der Anzahl der zugelassenen Zwischenmoden und 
wird auch als Auflosung oder Abstufung der Interpolation 
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zwischen den reinen Moden bezeichnet. Die Auflosung R kann 
einen beliebigen Wert annehmen, wird jedoch anwendungsabhangig 
zur Erzielung einer optimalen Systembeschreibung (insbesondere 
bei stark verrauschten Vorgangen) und praktikabler 
Rechenzeiten insbesondere unter Beriicksichtigung der oben 
genannten Schaltrate ausreichend niedrig gewahlt. Bei 
praktischen Anwendungen (siehe unten) kann vorgesehen sein, 
dafi die Auflosung R manuell von einem Bediener oder 
automatisch von einem Stellkreis in Abhangigkeit von einem 
vorliegenden Analyseergebnis und einem Vergleich mit einem 
vorgegebenen Schwellwert gewahlt werden. 

Die Gesamtzahl der gemischten Moden betragt bei einer 
gegebenen Auflosung R zwischen jeweils zwei Netzen 
|M | = R * N * (N-l)/2. Bei dem o. a. Beispiel betragt somit bei 
N = 8 reinen Moden und einer Auflosung R = 32 und Gesamtzahl 
der gemischten Moden |M| = 896. Zur Bestimmung der Gesamtzahl 
der Systemmoden kommen noch die 8 reinen Moden hinzu. 

Die Driftsegmentierung umfafit nun die Suche nach einer 
Segmentierung mit den reinen und gemischten Systemmoden (a, b, 
R) , die in Bezug auf den Vorhersagef ehler der Moden der 
gesamten Zeitreihe optimiert ist. Die Prediktoren werden so 
ausgewahlt, dafi jedem Element der Zeitreihe einer der Moden 
aus der Gesamtzahl der Systemmoden zugeordnet werden kann. Der 
Vorhersagefehler ist die Abweichung einer 
Prediktorenvorhersage vom tatsachlichen Element der zu 
untersuchenden Zeitreihe. Fur die zu untersuchende Zeitreihe, 
die nicht mehr unbedingt die Trainings-Zeitreihe ist, mit der 
bei der Schaltsegmentierung die angepaliten Netzwerke oder 
Prediktoren ermittelt wurden, wird fur jeden Zeitschritt mit 
jedem der Prediktoren eine Vorhersage ermittelt, woraus sich 
eine zeitabhangige Matrix der Prediktorerfvorhersagen ergibt, 
aus der ein mittlerer Vorhersagefehler fur willkurlich 
gewahlte Segmentierungen ableitbar ist. Die Segmentierungen 
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mit dem geringsten Vorhersagef ehler ist die gesuchte 
Drif tsegmentierung. 

Die Suche nach der Segmentierung mit dem geringsten 
Vorhersagef ehler kann mit jeder geeigneten Such- oder 
Iterationstechnik erfolgen. Vorzugsweise wird eine dynamische 
Programmiertechnik gewahlt, die Equivalent zum 

Viterbi-Algorithmus fiir HM-Modelle (sog. Hidden Markov Models) 
ist. Einzelheiten hierzu sind beispielsweise in der 
Publikation "A Tutorial on Hidden Markov Models and Selected 
Applications in Speech 'Recognition" von L. R. Rabiner in 
"Readings in Speech Recognition" (Hrsg. A. Waibel et al., San 
Mateo, Morgan Kaufmann, 1990, S. 267-296) beschrieben. Die 
Driftsegmentierung ist im Rahmen der HM-Modelle die 
wahrscheinlichste Modensequenz, die die zu untersuchende 
Zeitreihe erzeugt haben konnte. Als Nebenbedingung wird dabei 
die Moglichkeit der Moden-Anderungen iiber die T-Funktion 
(siehe unten) eingeschrankt . 

Das Ziel der Anpassung ist die Angabe einer optimalen Sequenz 
von Netzwerken bzw. linearen Mischungen aus diesen. Eine 
Sequenz ist optimal, wenn die sogenannte Energie- oder 
Cost-Funktion C* der Vorhersage minimiert ist. Die 
Cost-Funktion C* setzt sich aus der Summe der quadratischen 
Fehler der Vorhersage und den Cost-Funktionen der 
Modenubergange der Sequenz zusammen. Die Ableitung der 
Cost-Funktion C* jeweils zwischen zwei Zeitpunkten t 0 und t^ 
erfolgt induktiv, indem zunachst von einer Start-Cost-Funktion 
gemafi Gleichung (4) ausgegangen wird. 

C 3 (t 0 ) = e 3 (t 0 ) 

wobei 

e 3 (t) = (x t - gjlt^)) 2 

der quadratischen Fehler der Vorhersage der reinen oder 
gemischten Moden g ist. 



(4) 
(5) 
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Ftir den Induktionsschritt von t - 1 auf t wird die 
Cost-Funktion gemafi Gleichung (6) fUr alle s e S berechnet: 

C s (t) = e 3 (t) + min {C s (t-1) + T(s, s)}, t=t 0 +l, . . . tmax (10) 

Dabei ist T(s, s) die Cost-Funktion des Ubergangs von einer 
Mode s zu einer Mode s. 

Die optimale (minimale) Cost-Funktion C* ist dann: 
C* = min {C s (t max )} 

Die Funktion T entspricht bei den HM-Modellen den 
Obergangswahrscheinlichkeiten und kann anwendungsabhangig 
geeignet gewahlt werden. Es kann beispielsweise vorgesehen 
sein, abrupte Schaltubergange und gleitendes Driften zwischen 
zwei Netzen zuzulassen und alle anderen ttbergange durch T = oo 
auszuschliefien. 

Die Driftsegmentierung ergibt sich aus der ermittelten 
optimalen Sequenz von Netzwerken bzw. linearen Mischungen aus 
diesen, indent die Moden, die C* ergeben, zuriickverfolgt und in 
Abhangigkeit von der Zeit erfafit werden. 

Im AnschluJi an die Driftsegmentierung kann sich ein 
zusatzlicher Schritt der Reduzierung der Zahl der zur 
Modellierung verwendeten Netzwerke anschliefien, der im 
einzelnen unten beschrieben wird. 

Abschliefiend findet eine Identif izierung der segmentierten 
Moden statt, indem jedem Prediktor oder Vorhersagemodell die 
zugehbrige Systemmode zugeordnet wird. Diese Identif izierung 
erfolgt anwendungsabhangig. 
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Das Ergebnis der Drif tsegmentierung im Fall der oben unter 
Bezug auf Figur 1 erlauterten chaotischen Zeitreihe {x t } mit 
vier Moden wird im folgenden unter Bezug auf Figur 2 
beschrieben. Die Drif tsegmentierung umfafit die Suche nach 
einem Verlauf a(t), der einen speziellen Pfad zwischen den 
reinen Moden liefert, fur den der Vorhersagefehler der 
gesamten Zeitreihe optimiert ist. 

Jeweils nach den ersten 50 Zeitschritte mit der Mode gemafl f, 
erfolgt fur 50 Schritte ein zeitlinearer Obergang in die Mode 
gemaJi f 2 . Der Obergang ist ein zeitabhangiges Drif ten gemaB den 
Gleichungen (12) : 

f("x t ) = (1 - a(t))f^lt t ) + a(t)f 2 (lf t ) 

mit (12) 
a(t) = t b - ta fc a = 50, t b = 100 

Entsprechende Obergange finden nach dem 150., 250. und 350. 
Schritt fur jeweils 50 Schritte statt. 



In Figur 2 ist die Besetzung der jeweiligen Moden entsprechend 
den ermittelten Netzwerken in Abhangigkeit von der Zeit 
(Zeitschritte [1200, 2400]) dargestellt. Aus 
Obersichtlichkeitsgrunden sind die Obergangs- oder 
Driftbereiche entsprechend ihren Zeitgrenzen und Ausgangs- 
bzw. Endmoden in Rahmen dargestellt, in denen der jeweilige 
Driftverlauf zwischen den Moden gepunktet ist. Figur 2a zeigt 
bei einer Auflosung R = 32 (siehe Gleichung (3)) Obergange wie 
z. B. bei den Zeitschritten 1350 bis 1400 zwischen den Netzen 
2 und 4. Die Obergange sind wie gemafl Gleichung (12) erwartet 
linear. Bei einer geringeren Auflosung von R = 3 ergibt sich 
die in Figur 2b da.rgestellte Segmentierung. Abweichend vom 
linearen Driftverlauf sind die gepunkteten Obergange 
stufenformig. Dennoch ist die Darstellung auch bei der 
geringeren Auflosung eine adequate Beschreibung des 
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dynamischen Verhaltens des Systems, wie der Vergleich der 
zeitlichen Lage der Moden und des Driftens zeigt. 



(2) Anwendungsbeispiele fur die Erfassung von Dirif tilbergangen 

(i) Blutzellenregulation im menschlichen Kbrper 

Die Blutzellenregulation im menschlichen Korper stellt ein 
hochdimensionales chaotisches System dar, das durch die 
Mackey-Glass-Verzogerungs-Differentialgleichung (13) (siehe 
auch in der o. a. Abhandlung von J. Hertz et al.) beschrieben 
werden kann. 

dx(t)/dt = -0.1x(t) + (0.2x(t-t d ) / (l+x(t-t d ) 10 ) (13) 

Erfindungsgemafi konnen Zeitreihen von physiologischen 
Parameters die fur die Menge der roten Blutkorper 
charakteristisch sind, anwendungsabhangig segmentiert werden. 
Die Funktionsfahigkeit der Segmentierung wird im folgenden 
beispielhaft erlautert . 

Bei Vorgabe von zwei Moden A und B, die sich jeweils durch 
verschiedene Verzogerungsparameter t d = 17 bzw. t d = 23 
unterscheiden, erfolgt mit einer Abtast-Zeitschrittgrofie x = 6 
nach 100 Schritten zunachst ein Obergang von A nach B. Der 
Obergang dauert 100 Schritte und ist eine Oberlagerung der 
Gleichung (13) mit den beiden Verzogerungsparameter t d wahrend 
der Integration von Gleichung (13). Die Oberlagerung wird mit 
einem exponentiellen Drif tparameter a (siehe Gleichung (2)) 
gemafi Gleichung (14) erzeugt. 



a(t) = exp(-4t/100) , t = 1,..., 100 



(14] 



In der Folge wiederholen sich alle 100 Schritte stationare 
Moden A oder B bzw. die jeweiligen Obergange. FUr jeden 
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Riick-Obergang nach einem Drif tiibergang wird ein schaltartiges 
Umspringen angenommen. Figur 3a zeigt die entsprechende 
Zeitreihe iiber 300 Schritte. Die Driftsegmentierung mit sechs 
Prediktoren auf der Grundlage von RBF-Netzwerken mit jeweils 
40 Basisfunktionen, einem Einbettungsparameter m = 6 und dem 
Verzogerungsparameter i = 1 (siehe Gleichung (2)) ergibt das 
in Figur 3b gezeigte Bild. Es zeigt sich die erwartete 
Segmentierung der Zeitreihe in stationare Moden und 
Driftiibergange. 

Allerdings zeigt sich, dafl sich jeweils 2 Netzwerke auf eine 
Mode spezialisiert haben (2, 3 => Mode A, 5, 6 => Mode B) . In 
einer solchen Situation kann es erf indungsgemafi vorgesehen 
sein, den zusatzlichen Schritt der Reduzierung der Zahl der 
zur Modellierung verwendeten Netzwerke vorzusehen. 

Der Reduzierungschritt umfafit eine sequentielle Reduzierung 
der Netzwerkzahl, jeweils verbunden mit einer Bestimmung des 
mittleren Vorhersagef ehlers . Die Reduzierung (Abzug von 
Redundanznetzwerken) wird beendet, falls eine weitere 
Verringerung der Netzwerkzahl eine signifikante Erhohung des 
Vorhersagef ehlers bedeutet. Figur 3c zeigt das Ergebnis einer 
solchen Reduzierung. Der mittlere • quadratische Fehler RMSE 
bleibt bei der Reduzierung urn ein, zwei, drei und vier 
Netzwerke jeweils konstant, wohingegen ein starker Anstieg fur 
den Fall auftritt, dafi nur noch mit einem Netzwerk modelliert 
wird. Dies bedeutet, dafi das System optimal mit einer Zahl von 
Netzwerken modelliert wird, die gleich der Gesamtzahl von 
betrachteten Netzwerken, vermindert urn die Zahl von 
Redundanznetzwerken ist. 



Die adequaten Modellnetze werden dann erhalten, indem der 
RMSE-Wert fur jede Netz-Kombination mit reduzierter Netzzahl 
berechnet wird. Die Netz-Kombination mit dem geringsten 
RMSE-Wert umfafit die gesuchten Modellnetzwerke oder 
Prediktoren. Figur 3d zeigt die Driftsegmentierung nach dem 
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Reduzierungsschritt. Demnach beschreiben die ttbrig gebliebenen 
Prediktoren 2 und 5 das System vollstandig. 

(ii) Schlafdatenerfassung 

Eine weitere Anwendung der Erfindung liegt im Bereich der 
Analyse physiologischer Daten, die fttr den Ablauf von Schlaf- 
und Wachmoden von Lebewesen charakteristisch sind. Als 
Grundlage fttr sich anschliefiende Prozeduren zur Erfassung von 
Schlaf storungen konnen Zeitreihen z. B. von EEG-Daten 
segment iert werden. 

Figur 4a zeigt im Vergleich die Ergebnisse einer herkommlichen 
Schaltsegmentierung (oben) , einer Drif tsegmentierung (Mitte) 
und einer "manuellen" Segmentierung (unten) eines 
medizinischen Fachmannes (Schlaf forscher) auf der Grundlage 
von Erfahrungswerten am Beispiel eines Nachmittagsschlaf es 
eines gesunden Menschen. Die Schalt- und Drif tsegmentierungen 
erfolgen mit acht Netzwerken (netl ... net8) an 
Einkanal-EEG-Daten x(t) (Figur 4b). In Figur 4a sind wie bei 
Figur 2 aus Obersichtlichkeitsgrunden Rahmen gezogen, die bei 
den Driftmoden verdeutlichen, zwischen welchen Netzen 
interpoliert wird. Die gepunktete Linie im Inner en der Rahmen 
zeigt den jeweils tatsachlichen Verlauf . Die manuelle 
Segmentierung basiert auf der Beobachtung von physiologischen 
Signalen (z. B. EEG, EOG, ECG, Puis, Blutdruck, Atmung, 
Augenbewegungen) . Die Moden Wl, W2 bezeichnen zwei Wachmoden 
mit offenen bzw. geschlossenen Augen und die Moden SI, S2 
jeweils Schlaf zustande . "n.a." und "art." beziehen sich auf 
nicht in Betracht gezogene Zustande bzw. Artefakte. 

Die Schaltsegmentierung zeigt ein vergleichsweise 
undif ferenziertes Bild, das nur grob mit den Ubrigen 
Beobachtungen konsistent ist. So tritt beispielsweise ein 
Einschlafvorgang in alien drei Fallen bei t « 7000 auf. Die 
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Driftsegmentierung ergibt jedoch mehrere DriftUbergange, die 
zusatzliche Einzelheiten des Schlafverhaltens darstellen. Der 
"manuell" beobachtete Schlafbeginn zur Zeit t « 4000 wird 
durch einen exponentiellen Drif tUbergang von Netz net7 
(Wachmoden-Prediktor) zum Netz net4 (Schlafmoden-Prediktor) 
reprasentiert. Das Aufwachen setzt bei t * 9000 durch ein 
leichtes RUckdriften zum Netz net7 ein, das bis zum Erreichen 
des "manuell" ermittelten Aufwachpunktes t » 9500 gehalten 
wird. In dieser Situation erfolgt eine plotzliche Anderung des 
Wichtungsfaktors, so daii das Netz net7 eine groJiere Wichtung 
erhalt. Nach t « 9800 (Augen of fen) besteht eine Mischung der 
beiden Wachmoden-Prediktoren net7 und net2. 

(iii) Weitere Anwendungen und Vorteile 

Figur 4a zeigt, daft mit dem erf indungsgemaflen Verfahren 
automatisch detaillierte Segmentierungen erzielt werden 
konnen, die bislang nur durch die Beobachtung komplexer 
Merkmalsbilder auf der Grundlage breiter Erfahrungen und 
Intuitionen zuganglich waren. Dieser Vorteil ist nicht nur in 
der Medizin, sondern auch in anderen Gebieten nutzbar, in 
denen grofie Datenmengen bei der Beschreibung komplexer 
dynamischer Systeme anf alien. Derartige Gebiete sind die 
physikalische, chemische und/oder biologische 
Verfahrenstechnik, die Geologie, Meteorologie, Klimatologie, 
die Spracherfassung u. dgl . . 

Erfindungsgemalie Verfahren haben die folgenden Vorteile. Das 
betrachtete System kann hochdimensional sein (10 oder mehr 
Dimensionen) . Die Erfindung erlaubt eine Reduzierung der 
Komplexitat eines solchen Systems durch Betrachtung von 
niedrigerdimensionalen Moden und von wechselnden Obergangen 
zwischen diesen. Der Einsatz von Vorhersagemodellen fiir die 
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Segmentierung 1st invariant gegen Anderungen der Amplitude 
erfaiiter Signale. 

Die Anwendung der Erfindung fur die Vorhersage oder Kontrolle 
eines Systems erfolgt derart, daB zunachst wie oben 
beschrieben aus der vergangenen Beobachtung und der Kenntnis 
der aktuellen Moden der IST-Zustand des Systems erfaBt wird, 
der ggf. eine Mischung entsprechend dem Ergebnis der 
Driftsegmentierung darstellt. Der IST-Zustand entspricht einem 
dynamischen System f. Die Vorhersage bedeutet, daB das System 
f auf den momentanen Zustand x angewendet wird und sich daraus 
die Vorhersage fur den unmittelbar folgenden Zustand y ergibt. 
Die Kontrolle bedeutet, daB aus dem IST-Zustand die Abweichung 
von einem SOLL-Zustand ermittelt und aus der Abweichung eine 
geeignete Regelstrategie abgeleitet wird. 

Der Vorteil der Vorhersage und Kontrolle besteht darin, dafi in 
komplexen Systemen (z. B. bei der Erfassung chemischer 
Reaktionen in einem Reaktor) , die ggf. die Messung lediglich 
weniger MeBgrofien erlauben, die fur sich allein aufgrund von 
Mehrdeutigkeiten oder wegen systemimmanenten Verzogerungen 
keine direkten Ruckschlusse auf den Systemzustand und ggf. 
bestehende gemischte Zustande erlauben, dennoch detaillierte 
Informationen uber das System abgeleitet werden konnen. So 
kann bei dem Beispiel mit einer chemischen Reaktion aus der 
erfindungsgemaBen Erfassung bespielsweise der makroskopischer 
thermodynamischer ZustandsgroBen eine optimale Regelstrategie 
abgeleitet werden, die die Dosierung bestimmter 
Reaktionspartner umfafit. 
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Patentansprilche 

1. Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden Si , die jeweils einen Satz 
a(t) charakteristischer Systemparameter besitzen, wobei eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellierung unterzogen wird, die dazu eingerichtet ist, in 
jedem Zeitabschnitt einer vorbestimmten Mindestlange fur jede 
Systemvariable x(t) ein vorbestimmtes Vorhersagemodell f t fur 
eine entsprechende Systemmode Si zu erfassen, 
dadurch gekennzeichnet , dafl 

nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dem das System von 
einer ersten Systemmode Si zu einer zweiten Systemmode s 3 
tibergeht, eine Folge von gemischten Vorhersagemodellen g t 
erfaiit wird, die durch eine lineare, paarweise Oberlagerung 
der Vorhersagemodelle f t<j der zwei Systemmoden s i(j gegeben ist. 

2. Verfahren gema6 Anspruch 1, bei dem die Modellierung eine 
Schaltsegmentierung ist. 

3. Verfahren gemali Anspruch 2, bei dem die 
Schaltsegmentierung durch Simulierung einer 
Trainings-Zeitreihe des Systems oder der zu untersuchenden 
Zeitreihe mit mehreren, untereinander in Konkurrenz stehenden 
Vorhersagemodellen erfolgt. 

4. Verfahren gemafl Anspruch 3, bei dem die Vorhersagemodelle 
durch neuronale-Netzwerke oder andere Modelle zur Schatzung 
von Funktionen gebildet werden, die jeweils far eine Mode s 
charakteristisch sind und nach vorbestimmten Trainingsregeln 
urn die Beschreibung der einzelnen Elemente der Zeitreihe 
konkurrieren. 
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5 Verfahren gemafi einem der Anspriiche 1 bis 4, bei dem die 
Folge gemischter Systemmoden g, aus den Vorhersagemodelle f, . 
und Interpolationsparametern a, b gemafi 
g± = a(s)f i(s) (x) + b(s)f j(s) (x) ermittelt wird. 

6. Verfahren gemafi Anspruch 5, bei dem die 

Interpolationsparameter gemafi 0 < a(s) < 1 und b(s) = 1 - a (s) 
gewahlt sind. 

7. Verfahren gemafi Anspruch 6, bei dem die Werte a(s) auf 
eine bestimmte Auf losungszahl R beschrankt und/oder 
aquidistant sind. 

8. Verfahren gemafi einem der vorhergehenden Anspriiche, bei 
dem die Folge gemischter Vorhersagemodelle g t erfafit wird, 
indem fur jeden Zeitschritt mit jedem der moglichen 
Vorhersagemodelle jeweils eine Vorhersage ermittelt wird, 
woraus sich eine zeitabhangige Vorhersagematrix ergibt, aus 
der ein mittlerer Vorhersagef ehler fur willktirlich gewahlte 
Segmentierungen ableitbar ist, wobei die gesuchte Folge 
gemischter Vorhersagemodelle g t die Segmentierung mit dem 
geringsten Vorhersagef ehler bzw. der maximalen 
Wahrscheinlichkeit ist. 



9. Verfahren gemafi Anspruch 8, bei dem die Suche nach der 
Segmentierung mit dem geringsten Vorhersagef ehler mit einer 
dynamische Programmiertechnik erfolgt, die aquivalent zum 
Viterbi-Algorithmus fur Hidden-Markov-Model le ist, wobei eine 
optimale Sequenz von Vorhersagemodellen unter Verwendung einer 
minimierten Cost-Funktion C* der Vorhersage ermittelt wird und 
die Segmentierung induktiv aus der Sequenz von 
Vorhersagemodelle abgeleitet wird. 

10. Verfahren gemafi einem der vorhergehenden Anspriiche, bei 
dem nach der Drif tsegmentierung ein zusatzlicher Schritt der 
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Reduzierung der Zahl der zur Modellierung verwendeten 
Vorhersagemodelle erfolgt, bei dem die Zahl der 
Vorhersagemodelle sequentiell jeweils verbunden mit einer 
Bestimmung des mittleren Vorhersagef ehlers so weit reduziert 
wird, daB eine weitere Verringerung der Zahl der 
Vorhersagemodelle eine Erhohung des Vorhersagef ehlers 
bedeutet . 

11. Verfahren gemaii einem der vorhergehenden Anspruche, bei 
dem die Zeitreihe mindestens einer der Systemvariablen x(t) 
eine Zeitreihe von physiologischen Parametern umfafit, die mit 
der Mackey-Glass-Verzogerungs-Dif f erentialgleichung 

dx(t)/dt = -O.lx(t) + (0.2x(t-t d )/(l+x(t-t d ) 10 ) 
beschrieben wird. 

12. Verfahren gemaJi einem der Anspruche 1 bis 11, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von physiologischen Parametern umfafit, die ftir den 
Ablauf von Schlaf- und Wachmoden charakteristisch sind. 

13. Verfahren gemafl Anspruch 12, bei dem die physiologischen 
Parameter EEG-Signale umfassen. 

14. Verfahren gemaii einem der Anspruche 1 bis 10, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von Sprachsignalen umfafit. 
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Zusammenf as sung 



Bei einem Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden, die jeweils einen Satz 
a(t) charakteristischer Systemparameter besitzen, wird eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellierung wie z. B. einer Schaltsegmentierung unterzogen, 
die dazu eingerichtet ist, in jedem Zeitabschnitt einer 
vorbestimmten MindestlSnge fiir jede Systemvariable x(t) ein 
vorbestimmtes Vorhersagemodell wie z. B. ein neuronales 
Netzwerk fiir eine entsprechende Systemmode zu erfassen, wobei 
nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dent das System von 
einer ersten Systemmode zu einer zweiten Systemmode iibergeht, 
eine Folge von gemischten Vorhersagemodellen erfaflt wird, die 
durch eine lineare, paarweise Oberlagerung der 
Vorhersagemodelle der zwei Systemmoden gegeben ist. 
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